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Abstract. The interaction of a quantum field with a background containing a Dirac delta 
function with support on a surface of codimension 1 represents a particular kind of matching 
conditions on that surface for the field. In this article we show that the worldline formalism can 
be applied to this model. We obtain the asymptotic expansion of the heat-kernel corresponding 
to a scalar field on in the presence of an arbitrary regular potential and subject to this 

kind of matching conditions on a fiat surface. We also consider two such surfaces and compute 
their Casimir attraction due to the vacuum fiuctuations of a massive scalar field weakly coupled 
to the corresponding Dirac deltas. 

FAGS numbers: ll.lO.-z 

1. Introduction 

The Worldline Formalism (WF) is an intuitive and efficient method to compute effective 
actions in Quantum Field Theory (QFT) [1, 2]. This formalism has been applied to the coupling 
of external gravity to quantum fields of spin [3], spin 1/2 [4], spin 1 and, more generally, 
antisymmetric tensor fields [5]. Higher spin fields have also been studied in the WF for flat 
space [6] and a conformally flat background [7]. The method has also proved useful for the 
study at the quantum level of photon-graviton mixing in an electromagnetic background [8] . 

In the last years, the first steps of a systematic generalization of the WF to manifolds with 
boundaries have been taken through the computation of the asymptotic expansion of the heat- 
kernel trace for a scalar field on a flat manifold with boundary [9, 10, 11, 12]. 

Worldline technics combined with Monte Carlo methods have also been successfully applied 
to compute Casimir energies for various geometries [13, 14, 15], giving results on the edge effects 
[16] as well as an interesting behavior of the Casimir force in terms of the temperature [17, 18, 19]. 
One of the advantages of worldline Monte Carlo with respect to numerical calculation on the 
lattice is that the former can be implemented on a fermionic system without a discretization of 
spacetime [20]. 

In this article we show that the WF can be also applied to a field under the matching conditions 
on a surface of codimension 1 imposed by a background containing a Dirac delta with support on 
that surface. To do that, we obtain the asymptotic expansion of the heat-kernel corresponding 
to a scalar field on 1^"^+^ interacting with a delta function with support on a flat d-dimensional 
surface. Later, in the same formulation, we study the Casimir force between two such surfaces 
due to the vacuum oscillations of the scalar field. 

A Dirac delta potential with support on a surface of codimension 1 is a well-defined mathemat- 
ical problem [21] which has been extensively studied to model several physical configurations. 

In the last years, this setting has been applied to quantize the TE-modes of the electromagnetic 
field in the vicinity of systems such as giant carbon molecules [22, 23, 24, 25, 26, 27, 28, 29, 30]. 
The collective electrons of these molecules can be considered as an infinitely thin plasma layer 
that imposes semi-transparent boundary conditions on the oscillation modes of the quantum 
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electromagnetic field. For small values of the net charge and current at the layer the problem 
reduces -for TE-modes- to solving the wave equation with a delta-like potential at the layer. 

Effective actions for quantum fields at the 1-loop level can be computed using zeta-function 
and heat-kernel technics [31, 32]. In particular, the asymptotic expansion of the heat-kernel for 
small values of the proper time contain information of the 1-loop ultraviolet divergencies. In 
the WF this asymptotic expansion corresponds to the semiclassical expansion of a field theory 
in 0+1 dimensions. 

Let us consider a scalar field 4>{^) defined on a: = {xi, . . . , Xd+i) G with Euclidean action 
(1) S[ct>]= [ \ld,<P{x)d,c^ix) + U{cP{x)) 



2 

The corresponding 1-loop effective action is given by 

(2) m = S[4>] + hogBetA, 

where the operator of quantum fluctuations A -for the action (1)- is the differential operator 

(3) A = -A + V{x) := -A + U"{(t>{x)) . 

The primes on U {(p) denote derivatives with respect to (j)- The functional determinant of the 
Schrodinger operator A can be defined by means of its corresponding heat-kernel K. A usual 
definition for the functional determinant in expression (2) is [32] 

(4) log BetA:= - r ^ Tr e-^^ , 

Jo P 

where 

(5) Tre"^^ = / K{x,x,^) 

and the heat-kernel K is defined as the solution of the system 

(6) {di3 + A)K{x,x',p) = 0, K{x,x',0) = 6{x,x') . 

The divergencies at /3 = of the integral in eq. (4) are related to the ultraviolet divergencies 
of the effective action. Consequently, the regularization of these divergencies can be carried out 
by performing the integral in /3 from to oo, where A plays the role of an ultraviolet cut-off 
regulator. 

For a differential operator of the type (3) defined on a compact {d + l)-dimensional base 
manifold, it is a well-known result of Spectral Geometry [33] that -under quite general hypothesis 
on the manifold, the boundary conditions and the regularity of the potential- the heat-kernel 
trace Tr e~^^ admits the following asymptotic expansion for small values of the proper time /3 

^ oo 

where the Seeley-de Witt coefficients 6„/2 integrals on the base manifold and its boundary of 
local invariants built from the potential V and the geometry of the base manifold. If, as already 
mentioned, one uses a cut-off regularization of the integral in eq. (4) then, as can be seen from 
eq. (2), the first coefficients 6„/2 give the 1-loop counterterms that regularize the effective action. 
It is important to mention that the coefficients 6„/2 foi' odd n, i.e. those multiplying half-integer 
powers of /3, vanish if the base manifold has no boundary. 

A different asymptotic expansion of the heat-kernel could be obtained if one relaxes some of 
the regularity hypothesis required for the validity of (7); many examples along this line have 
been extensively studied (see the review [32]). In particular, the asymptotic expansion of the 
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heat-kernel trace for a Schrodinger operator whose potential contains a Dirac delta function has 
been determined, up to our knowledge, for the first time in [34]. 

The purpose of this article is to show that the worldline formalism can be used as a very 
efficient tool to study this kind of singular potentials. In Section 2 we will review some basics 
of the WF by computing the asymptotic expansion of the heat-kernel K{x,x',/3) for small /3 
corresponding to a Schrodinger operator A given by (3) with an arbitrary regular potential V{x). 
We will then replace this expansion into expression (5) to obtain, according to expansion (7), 
the first Seeley-de Witt coefficients. In Section 3 we will generalize the WF for the case in which 
the potential is not regular but contains a Dirac delta with support on a flat layer. Finally, in 
Section 4 we use these techniques to compute the Casimir energy between two such layers in the 
weak coupling limit. 



2. Regular potential 

In this section we will show how the WF works by computing the first terms in the asymptotic 
expansion of the heat-kernel corresponding to the Schrodinger operator (3). Moreover, in the 
next section this result will prove useful to extend the WF to the case where the potential 
contains terms with Dirac deltas with support on flat surfaces. 

First of all, notice that according to eqs. (6), the heat-kernel K{x, x' , /3) is the wave function at 
X of a particle that was initially at x' and has evolved an Euclidean time f3 with Hamiltonian A = 
—A+V{x) . This is the first step in the implementation of the WF: since the operator of quantum 
fiuctuations A that determines the 1-loop effective action (see eq. (2)) depends on the field theory 
under consideration, one must find for each field theory an auxiliary quantum mechanical particle 
with Hamiltonian A, such that the corresponding heat-kernel could be regarded as its wave 
function. 

In the case of the self-interacting scalar field with action (1) the operator A is given by (3) 
and therefore corresponds to the Hamiltonian of a non-relativistic Schrodinger particle. The 
next step in the WF consists in implementing the path integral quantization of this auxiliary 
particle: 

(8) K{x, x\ fi)= j Vx{r) e- [i\t)/4 + V {x{t))] 

x{0)=x' 
x{l3)=x 

In order to get an expansion for small /?, it is convenient to make the following rescaling in 
Euclidean time r — )• /3r so that 



(9) 



i^(x,x',/3)= / Vx{r)e-!idt[i\t)/AP + pV{x{t))] 



x(Q)=x' 
x{l)=x 



This expression is indeed appropriate for a small /3-expansion for two reasons. On the one hand, 
since j3 acts as the intensity of the potential term, one gets a series in powers of V 



(10) K{x,x',^) = ±tPL j p,(,),-/;dt'x^(W 



x{0)=x' 
x{l)=x 



1 

dtV{x{t)) 







On the other hand, the kinetic term shows that the propagator is proportional to /3, so a 
complete expansion for small /3 can be flnally obtained by a loop expansion in the sense of the 
-|- 1 quantum field x{t). In order to do that we consider the quantum field x[t) in terms of 



4 S.A. FRANCHINO VINAS AND P.A.G. PISANI 

the fluctuations <^{t) about the classical trajectory with respect to the case V = {) 

(11) Xi{T) = {Xi - x[) T + x[ + Lpi{T) . 

Notice that the perturbation fields '^{t) satisfy Dirichlet boundary conditions at r = and 
T = 1. Replacing these new fields (^(r) in eq. (10) straightforwardly gives 

rl\2 oo 



<^(1)=0 



/ dtV {{x-x')t + x' + ip{t)) 
Jo 



(12) =e 4/3 ^LJ^/I dtV{{x-x')t + x' + ip{t))] 

In the last equation, we have defined the mean value 

(13) (...):= / V^ir) e- fo dt' ^HtW ... , 



</p{0)=0 

v(i)=o 

where the dots represent a function of the fields (p{t). 

At this point, it is important to remark that expression (12) coincides with the free heat-kernel 
Kq{x,x' , 13) corresponding to V = 



(14) Koix,x',/3) 



(x-x'f 

e 4/3 



(4^^)(d+l)/2 ' 

if we normalize 

(15) (1):= / P^(r)e-/o^*V^(t')/4/3, 



■ (47r/3)('^+i)/2 ■ 

</p(0)=0 
¥'{l)=0 

The next step is to Taylor expand the regular potential V about x' to get an expansion in 
powers of ^fi{t) and the increments {xi — x'^. The first contributions coming from the terms 
n = 0, 1, 2 in the series of expression (12) read 

(16) K(x,2;',/3) = e 4/3 x 

X |(1) -/3y(rE')^' dt{l) - pdiV{x') dt [{xi-x'i)t{l) + {^i{t))] - 

-^didjV{x') f dt {[{xi-x'^)t + ^^{t)\[{xj-x'j)t + ^j{t)\)- 



...+^V\x')j^ / rftdt'(l) + .. j . 
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Each term of this expansion contains an n-point transition function 

. . . ^M) ■■ = I ^Vir) e- /o dt' ^'it')/il3 ^^^^^^^^ _ _ _ ^^^^^^^ 



</3(0)=0 



The generating functional Z[j] is defined as 

Z[j]:= I V^{r)e-fodt'mt')/^P+j{t'Mt')} 
<^(o)=o 

¥'(l)=0 

According to expressions (17) and (18), every n-point transition function can be computed in 
terms of the propagator, which is proportional to the one-dimensional symmetric Green function 
under the corresponding Dirichlet boundary conditions at f = and t = 1 

23 

(19) {^^m{t')) = 5., J^^^^^^ t{l - t') ift<t'. 

We can now compute the mean values in (16) and straightforwardly obtain the first terms of 
an expansion of the heat-kernel in powers of f3 and the increments 5i := Xi — x[ 

_|i , 



(20) 



T ^ df^^ . 5,5,6, + ^ (-^^AAF + VAV + \d,Vd,V - V^^ + 



where the upper and lower sign apply if the potential and its derivatives are evaluated at the 
initial point x' or final point x, respectively. Notice that, because of the Gaussian prefactor, the 
integration of the kernel K{x,x',l3) in x implies that each factor 6, would give a contribution of 
order /S^/^. Therefore, the dots in expression (20) represent terms of order /J''/^ . 

To obtain the asymptotic expansion of the heat-kernel trace we consider the heat-kernel at 
the diagonal x = x', given by expression (20) for 6, = 



P'' f 1 .... , ..... , I0..0.. ..3\ , ^,«4^ 



+^ |^-_AAy + yAy + -aiF5,F-F^j +o(/3^)j , 

where the potential and its derivatives are now evaluated at x. According to eqs. (5) and (7), 
we integrate expansion (21) in x € M'^"''^ to obtain the first Seeley-de Witt coefficients (see e.g. 
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[32]) 
(22) 



bo= [ 



h 



1 



1 

6 Jwd+l 

1 

60 



V{x) 

AV{x) + 3V^{x) 

-AAV{x) + 10 V{x) AV{x) + 5 d^V{x) diV{x) - 10 V^{x) 



Notice that, as already mentioned, since the base manifold M.'^^^ has no boundary, the asymptotic 
expansion of the heat-kernel trace contains only integer powers of f3. 

In the next section we will show how this procedure can be generalized to the case where 
the potential contains a Dirac delta term. In that case, a straightforward application of the 
method developed in this section would certainly break down at the point of eq. (16) for a 
Taylor expansion of the potential would lead to products of the Dirac delta and its derivatives 
evaluated at the same point. 



3. Delta function potential 

In this section we apply the WF to compute the asymptotic expansion of the heat-kernel 
trace corresponding to a Schrodinger operator whose potential contains a Dirac delta with 
support on the hyperplane xi = of IR'^+^, i.e. 

(23) A^ = -A + V{x)+j5{xi) , 

where x G ]R''+\ V{x) is a regular potential and 7 G M+. We will denote y G M'^ the coordinates 
on the surface xi = 0. 

The heat-kernel at the diagonal corresponding to the operator A^ is given by (see eq. (9)) 
K,ix,x,P)= [ Vx{r)e-Iodt'[iHt')/4^ + pV{x{t')) + P^6[x,it% ^ 



x{0)=x 
x{l)=x 

Notice that we have already rescaled r — >• /3r. In order to use the results of the preceding section 
we expand the exponential containing the delta funcion 

(24) K^{x,x,P) - Ko{x,x,(3) = ^ J ■ ■ ■ J dh • • ■ dtn x 



n=l 







ir)e-fodt' [i^(W + /3F(x(t'))]5[,^(,^)]...5[,^(,„)] . 



Vx 

x{0)=x 
x{l)=x 

In this expression we have subtracted the contribution Kq{x,x,I3) corresponding to the case 
7 = 0, which has been already computed in the previous section (eq. (21)). 

As already mentioned, now we can certainly not expand the delta functions in small fluc- 
tuations about the same point; otherwise, divergent products of delta functions would arise. 
Instead, we regard these path integrals as transition functions under the regular potential V{x) 
but constrained by the delta functions' terms. 
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Indeed, the nth. term in the series of expression (24) receives contributions from all closed 
"trajectories" at the point x with the restriction that at "times" ti, . . . t„ the "classical particle" 
should "hit" the surface xi = 0, which is the support of the delta function. Consequently, each 
path integral with n delta functions is computed as the product of n + 1 transition functions 
under the n constraints xi{ti) = . . . = xi{tn) = 0. Expression (24) can be therefore written as 

K^{x, x, /3) — Kq{x, x, (3) = 

!ldt'y\t')/AB 



Y.{-PlT j ■■■ j j dtidt2 ...dtn>^ j Vy{T)e 

"=1 y{Q)=y 

yW=y 

VxUr) e-^o^t' [x,\t')/Af3 + pV{x{t'))] 



, . xUUj=a:i 
(25) xi(ti)=0 



xi{ti)=0 



Cdt' [x,\t')/AP + pV{x{t'))\ 



j Vxiir) e 



a;i(tn)=0 
3;i(l)=a;i 

In this expression we can see that, for fixed trajectories ?/(r) in the coordinates parallel to the 
surface, the path integrals in xi{t) are related to one-dimensional heat-kernels in the presence of 
a regular potential. Notice, however, that the regular potential as a function of the variable xi 
depends explicitly on time, due to its dependence on the fixed trajectories ^(t). In spite of that, 
these heat-kernels could still be computed along the lines developed in Section 2. Afterwards, this 
procedure should be repeated on the path integrals on the trajectories 2/(r). Finally, integrating 
in a; £ W^~^^ one would obtain the asymptotic expansion of K{x,x,(3). 

Nevertheless, as a matter of technical simplicity, instead of using expression (25) we will 
make further use of this interpretation of the delta potentials in the path integral as constraints 
on the trajectories. It is more convenient to introduce for each delta function 6{xi{ti)) in 
the path integrals in (24), an extra d-dimensional delta function on the coordinates y, namely 

6^'^'^ {y{ti) — yQ^). These delta functions are then implemented as the constraint x{ti) = (0, y^^) as 

(i) 

already explained. Next, integration in the variables y^ renders the expected result. Expression 
(24) then reads 

(26) K^{x, X, /3) - Ko{x, x, (3) = 

1 h t2 

dy^^K . . dy^""^ X 



V(-/3 7)" [■■■[[ dhdt2 ...dtn X [ ... [ 

„_i J J J JR'* Jw 





X K{x^^\x,l3h) X ... i^(x(^+^),x«,/3(i,+i-i,)) ••• X i^(x,x("),/3-/3t„) , 

where x*-*^ = (0, y^*^). The asymptotic expansions of the kernels K in expression (26) can be 
read from expression (20). 

A remark is now in order. It could be argued that if the integrals in the coordinates y^*^ 
were absent and the coordinates j/^*-* took a fixed value, say y^*) = 0, then expression (26) would 
give the transition function in the presence of a d + 1-dimensional delta with support at the 
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origin. However, it is well-known that this kind of potential gives an ill-defined Schrodinger 
Hamiltonian. This is consistent with the fact that if the integrals in the y^*-* were absent in 
expression (26), then the subsequent integrals in ti would be divergent. 

Turning back to expression (26), after replacing the heat-kernels by their expansions given 
in (20) we can finally integrate in rr G M'^"'"^ to obtain the small /3-asymptotic expansion of the 
difference between the heat-kernel traces 

(27) Tre-^'^'-Tre-^-^ = ^^-^^ld%^^^ 

7 (8^ - 7^) /3^/' + + 7' V) (3' + 0(/3V2) | ^ 

where denotes the regular Schrodinger operator given by (23) for 7 = 0, whose heat-kernel 
expansion was computed in Section 2. 

Prom expansion (27) we read the corrections A&„ to the first Seeley-de Witt coefficients bn 
for n = 1,2,3 (see eqs. (22)) due to the Dirac delta superimposed to the regular potential V. 
Moreover, we see that the presence of the delta potential implies the appearance of half-integer 
powers of /3 in the heat-kernel expansion, as is the general case when the base manifold has 
boundaries. 

We display the corrections to the first Seeley-de Witt coefficients, as well as the first coefficients 
A6„/2) with odd n, corresponding to half-integer powers of /?: 



Abi 


= -7 




~ 4 ^ 


A62 




A65/2 


32 



(28) ^U2^ ryv -^7^ 

^■7' (8F-7') 

Ab3 = lj(AV -3V^ + ^^V) . 
6 

The first four coefficients displayed in (28) can be obtained as a particular case of the results in 
[34]. 

4. Casimir force between two semitransparent boundaries 

The calculations of Section 3 can also be carried out in the presence of Dirac delta functions 
with support on several hyperplanes. As mentioned in the introduction, a Dirac delta can be 
used to model the boundary conditions of the TE modes of the electromagnetic field in the 
vicinity of plasma layers for small values of the fluctuations of the charge density and current 
in the layer. Consequently, as an application, we will compute in this section the Casimir force 
between two parallel plasma layers due to the quantum oscillations of a scalar field which is 
weakly coupled to the Dirac deltas. This has already been computed for a general coupling by 
using an explicit expression of the corresponding Green function in [35]. 

We consider the following Euclidean action for a massive scalar field ipix) defined on x G M°'+^ 

(29) 5[0] = 1 [ di<P{x)di^{x) + m'^(P'^{x) + 7 [5{xi + L/2) + 6{xi - L/2)] (l)\x) , 



2 

where the interaction with the Dirac delta functions imposes semitransparent boundary con- 
ditions on the scalar field at the hyperplanes xi = —L/2 and xi = L/2. Notice that the 
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d + 1 coordinates in Euclidean space consist in the Euclidean time, the spatial coordinate xi, 
describing the position of the layers, and d — 1 spatial coordinates parallel to the layers. 
The corresponding heat-kernel at the diagonal is given by (see eq. (9)) 

(30) K{x, X, /3) = e-l^""' x 



j Vxiir) e 



dt' [xi^{t')/Af3 + /3 7 5 [xi (f) + L/2] + (3 j 6 [xi{t') - L/2]] 



xi{0)=xi 

Xl{l)=Xl 



(31) 



where we have factorized the contributions from Euclidean time and the d — 1 coordinates 
parallel to the layers. As we have mentioned, we are interested in the Casimir energy in the 
weak coupling regime, i.e. for small values of the coupling parameter 7. 

A rescaling of the dimensionful constants m and 7 in terms of the distance L between the 
layers leads us to the following expression for 7L <^ 1 

Kix,x,P) = ^^e-^-' I Vx,ir)e-fodt'^.HtW, 

xi{0)=xi 

Xlll) = Xl 

x|l-/37^ dt5[xi{t) + L/2]- p-i dtS[xi{t) - L/2] + 

+ / [ dsdt6[xi{s) + L/2]6[xi{t)-L/2] + . 
Jo Jo 

The first term in the series (31) is related to the divergent energy of empty space, while the 
second and third terms contribute to the self-energy of each layer, which is easily shown to be 
independent of the separation L between the layers. It is the fourth term that gives the leading 
contribution in 7 to the interaction energy between the layers. This term, which we denote by 
Kl{x,x, /3), can be computed as described in Section 3 and gives the following contribution to 
the heat-kernel trace 

dx Kl{x, X, f3) := 



TAd- 



(47r/3)'^/2 



1 ^2^2^ 



1 rl 



(32) 



dxi / Pxi(r) e-io'^*'^'i'(*')/'^^ / / dsdt6[xiis) + L/2] 5[xi{t)- L/2] 
J Jo Jo 

xi{0)=xi 

Xl{l) = Xl 

TA,., -^2^2 ^2^ 



(47r/3)^/2 

^ * g-(xi+L/2)2/4/3s g-LV4/3(t-^*) g-(xi-L/2)V4/3(l-t) 



/ 



dxi I dsdt , , , 

V4^F^ y/47r/3(i- s) v^47r/3(l -t) 





2 

erfc(L/v^) . 



(47r/3)^/2 4 

In this expression T represents the infinite length of the time interval, Aii_\ is the volume of the 
[d — l)-dimensional layers and erfc is the complementary error function. 
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The Casimir energy Eq per unit area of the layers can be defined by means of the 1-loop 
correction to the effective action r[i;^] (given by the second term in eq. (2)) as Eq := 
/T A(i~i. Therefore, the interaction energy per unit area between the layers due to the 



vacuum oscillations of the scalar field reads (see eq. (4)) 

(33) Eo = --^]- [ dxKLix,x,/3) , 

-t Ad^i 2 Jq P J^d+i 

whereas the Casimir pressure on the layers can be computed as 

Replacing the leading contribution of the interaction between the layers to the heat-kernel given 
by eq. (32) into eqs. (33) and (34) we obtain the pressure in the weak coupling regime <C 1 

(35) p = - ^^jl,y, {mlLfl^~'l^ K,/,^,/2{2mL) , 

in terms of the modified Bessel function. 

Let us finally mention two limiting cases of expression (35). Firstly, we consider the Casimir 
pressure for large separation L of the layers. If we take the limit mL 1 in eq. (35) we obtain 

2 „rf/2-l 

As remarked in [35], the exponential decrease with the product niL of the Casimir attraction 
due to the vacuum oscillation of massive fields justifies the consideration of photons for the 
determination of the Casimir effect. 

Secondly, we consider the massless limit of expression (34). For mL ^ 1 we obtain 

2 

J log {Lm) if d = 1 , 



(37) p ' 



f nid-l)/2) 1 

(4^)(d+i)/2 2 Ld~i II a ^ ^ • 



The first line of eq. (37) is consistent with the fact, pointed out in [36], that in 1 + 1 dimensions 
(d = 1) the massless field generates a Casimir pressure which is not analytic in 7 for small values 
of this parameter. For d = 3, the second line can be obtained as the weak coupling limit of the 
result in [26]. 

5. Conclusions 

In this article we have considered the generalization of the Worldline Formalism' (WF) tech- 
niques to study scalar fields under semi-transparent boundary conditions. In particular, we 
determine the asymptotic expansion for small values of the proper time of the heat-kernel corre- 
sponding to a Schrodinger operator whose potential contains Dirac delta functions with support 
on flat surfaces. 

In Section 3 we have studied the case in which the potential contains a regular term su- 
perimposed to a single delta function. In Section 4 we have considered two delta functions 
with support on parallel surfaces and we have computed the Casimir energy due to the vacuum 
oscillations of a scalar field weakly coupled to the delta functions. 

We have shown that the WF, which has been extensively used as a very efficient and intuitive 
tool for 1-loop calculations, is also very suitable for handling with delta functions. Indeed, some 
advantages of this formalism are still present in the case of a potential with delta functions. 
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On the one hand, although we have computed the small /3-asymptotic expansion of the heat- 
kernel trace, the WF provides a mechanism for computing the heat-kernel K{x,x',/3) out of 
the diagonal, i.e. for x 7^ x'; this determines, e.g., the expansion of the propagator at small 
distances. 

On the other hand, the techniques of the WF do not require the knowledge of the explicit 
solutions of the corresponding Schrodinger operator. For this reason, we were able to compute 
the effects of the delta function on the Seeley-de Witt coefficients corresponding to a regular 
potential. In relation to the aforementioned application of this model to the study of the Casimir 
energy of the TE modes in the vicinity of flat plasma layers, the calculation of the present article 
would be related to the case where the space between the layers were filled by a dielectric with 
position dependent-permittivity. A natural generalization of our result, would be to consider 
the case in which the Dirac delta functions have support on curved surfaces. 
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